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Abstract 

We demonstrate that necklace-shaped arrays of localized spatial beams can merge into stable 
fundamental or vortex solitons in a generic model of laser cavities, based on the two-dimensional 
complex Ginzburg-Landau equation with the cubic-quintic nonlinearity. The outcome of the fusion 
is controlled by the number of "beads" in the initial necklace, 2N, and its topological charge, M. 
We predict and confirm by systematic simulations that the vorticity of the emerging soliton is 
\N — M\. Threshold characteristics of the fusion are found and explained too. If the initial radius 
of the array (Ro) is too large, it simply keeps the necklace shape (if Rq is somewhat smaller, the 
necklace features a partial fusion), while, if Rq is too small, the array disappears. 
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I. INTRODUCTION 



The complex Ginzburg-Landau (CGL) equation is a universal model which plays an 
important role in many areas such as in superconductivity and superfluidity, fluid dynamics, 
react ion- diffusion phenomena, nonlinear optics, Bose-Einstein condensation, and quantum 
Se,d theories The CGL m ode, S produce dy— behaviors, e.g., spatiot_l 

chaos and formation of periodic patterns to dissipative solitons and their bound states 1, 

Many recent works focused on localized complex patterns in conservative models of op- 
tical media, different from the simplest nodeless ground-state modes, such as vortex soli- 
ions 6|, soliton clusters [3, lsj|, dipole-mode structures and their multipole counterparts 
si [lO, 11 , Q, Q], and necklace-ring solitons (NRSs) Q, [l5, Q. Similar patterns may 
be stable in dissipative models based on the CGL equation with the cubic-quintic (CQ) 
nonlinearity. These include stable localized vortices in two- and three-dimensional (2D and 
3D) CGL equation [l?], Q, Q and necklace-shaped soliton clusters 20 L Recently, stable 
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spatiotemporal NRSs were also reported in the 3D CQ-CGL equation 

Previous studies of complex dissipative-solitons patterns were focused on their stability. 
A more general issue is a possibility of dynamical transformation (fusion) of an initial soliton 
cluster into complex structures of a different type, or into a simple fundamental soliton. In 
this work, we demonstrate that necklace patterns can fuse into stable fundamental and vortex 
dissipative solitons in the 2D CQ complex GL equation, which is controlled by topological 
numbers of the initial pattern. Such an outcome is only possible in dissipative models, where 
the energy and angular momentum are not dynamical invariants. We predict the vorticity 
of the eventual state, and corroborate it by numerical results. If the initial radius of the 
necklace array is too large, it features a partial fusion, or simply keeps the initial structure; 
on the other hand, if the initial radius is too small, the pattern decays to zero. 



II. THE MODEL 

Following the notation adopted in Ref. [17J , we consider the CGL equation with the CQ 
nonlinearity in the (2+l)-dimensional setting, with propagation distance Z and transverse 
coordinates, X and Y: 
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iuz + idu + (1/2 — if3){uxx + uyy) + (1 — ie)\u\ 2 u — (i> — ifi)\u\ 4 u = 0, (1) 

where a > and /i > are parameters of the linear and quintic loss (the latter one 
accounts for the gain saturation), e > is the cubic gain, v (that may be both positive and 
negative) accounts for the quintic self-defocusing/focusing, the diffraction and cubic self- 
focusing coefficients are normalized to be 1, and (3 > is the effective diffusion coefficient. 
Except for the latter one, all other coefficients are standard ingredients of optical models 
based on the CGL equations; as for (3, it appears in models of laser cavities, as 
— t p t c A/(t p + r c ) 2 (in normalized units), where r p and r c are the polarization-dephasing and 
cavity-decay times, and A detuning between the cavity's and atomic frequencies. Therefore, 
the relevant case of (3 > corresponds to the negative detuning. 

We stress that the CQ nonlinearity in Eq. (1) is not a truncated expansion of a saturable 
nonlinearity, but represents a fundamental response of the nonlinear medium, due to some 
intrinsic resonance (s) in it. This response was directly observed in chalcogenide glasses 
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and m so n e orgamc opt,cal materials H - although saturable nonlinearity is more 

generic in optics [2j. 

The model based on the CGL equation (1) does not take into regard the finite relaxation 
time of the gain and loss in optical media, which may be added to the model through the 



respective evolution equation 
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271 ] . However, it seems quite plausible that the results 
reported below will not be strongly affected by this modification of the model. 

It has been demonstrated that vortex solitons and soliton clusters can self-trap in the 
framework of the 2D CQ-CGL model (and its 3D extension), provided that coefficient (3 is 



presented 0, Q, Q, ]20|, ]2l|. The stability of the localized patterns is supported by the 
simultaneous balance between the transverse diffraction and self-focusing, and between the 
gain and linear and quintic losses (including the effective diffusive loss, accounted for by 
P>0). 



Following Refs. jbll Q, [if]] , the initial necklace-shaped pattern, with amplitude A, mean 
radius Rq, and width w, can be taken (in polar coordinates r and 9) as 

u(Z = 0, r, 6) = Asech[(r - R )/w] cos(N6)exp{iM6), (2) 
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Here, integer N determines the number of elements ("beads") in the ring (necklace) struc- 
ture, which is 2N ', while another integer, M, is the topological charge of the complex pattern. 



III. NUMERICAL RESULTS 

Analysis of numerical results demonstrates that generic outcomes of the evolution of 
the necklace array can be displayed, e.g., for a = 0.5, (3 = 0.5, e = 2.5, v = 0.01, and 
jj, — 1, which corresponds to a physically realistic situation and, simultaneousl y m akes the 
evolution relatively fast, thus helping to elucidate its salient features 

Q, Q, Q- In this 

case, the amplitude and width of the individual 2D stable fundamental soliton, as found 
from Eq. (1), are A = 1.6 and w — 1. The robustness of the emerging patterns was tested in 
direct simulations of Eq. (1) with the initial condition taken as expression (2) additionally 
multiplied by [1 + p{x)} where p(x) is a Gaussian random function with zero average, the 
mean size of the perturbation amounting to 10 percent of the soliton's amplitude. 

The array with the topological charge equal to half the number of "beads" in the necklace, 
i.e., M = N, merges into a fundamental soliton, under the condition that initial radius Ro 
of the necklace array is smaller than a certain critical value, R^axi see Fig. 1(a). As 
follows from Eq. (2), the mean phase shift between adjacent beads in the initial array is 
Aip = Aip Q + 27rM/ (2N) = 0, where A(p Q = —n corresponds to the opposite sign of adjacent 
"beads" in the necklace with M — 0. Therefore, individual elements in the array, being 
in-phase, attract each other, which leads to their fusion into a stable fundamental soliton, 
as seen in Figs. l(b, c). 

Note that the CGL equation, being a dissipative one, does not have any dynamical 
invariant, hence the total angular momentum is not conserved either. In fact, the initial 
configuration corresponding to Eq. (2) with M = N may be realized as a mixture of two 
states, with values of the vorticity (total topological charge) S = M—N = and S = M+N. 
Obviously, the effective diffusion term in underlying Eq. (1), which is proportional to f3, 
produces a much stronger dissipation effect on the component of the wave field with the 
nonzero vorticity, hence only the zero- vorticity one survives in the limit of Z — > oo, as 
observed in Fig. 1. 

It is also easy to understand the increase of R^ ax with N, see Fig. 1(a). Indeed, the 
attraction between adjacent "beads", necessary for their fusion into the fundamental soliton, 
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is not too weak if the separation between them does not exceed a certain (maximum) value 

n 

[5[. The radius of the necklace, corresponding to a given separation between the "beads", 
grows linearly with N, which explains the roughly linear form of dependence _R^ ax (iV)in 
Fig. 1(a). 

When radius Rq of the initial necklace pattern is larger than a certain minimum value, 
R min , but smaller than 1.8N, i.e., 

R V mm <Ro< l-8iV, (3) 

the array with topological charge M evolves into a stable vortex soliton, provided that M 
falls into either of two mutually symmetric intervals, M^ in < M < N or N < M < M^ ax < 
2N — M^in see Fig. 2. As said above, initial configuration (2) may be realized as a mixture of 
states with vorticities whose absolute values are S = \N±M\, hence the vorticity component 
which may survive in the course of the evolution (one which is least affected by the diffusion 
term) is Sfi n = \N — M\. Indeed, the simulations confirm that the emerging vortex soliton 
features precisely this value of the vorticity. Because the asymptotic form of the vortex at 
r — > is const -r 3 , smaller values of S = \N — M\ correspond to a smaller radius of the 
inner hole in the vortex soliton, as observed in Figs. 2(c-e). Dependence M^ in (Ro), which is 
displayed in Fig. 2(b) for fixed values of N, can be explained too. Indeed, well-pronounced 
minima in the dependence at intermediate values of Ro corresponds to the fact that the 
fusion of the necklace array into vortex soliton is most feasible when the initial radius of the 
pattern is close to the radius of the expected vortex ring, which determines the location of 
the minima in Fig. 2(b). 

On the other hand, when the initial radius of the necklace is smaller than the minimum 
value, R min [see Eq. (3)], the pattern rapidly disappears (decays to zero) upon the propaga- 
tion, irrespective of the value of M , as shown in Fig. 2(f). This effect is explained by the 
strong dissipation generated by the effective diffusion term in Eq. (1) (the one proportional 
to 0) in the necklace array of a small radius. Naturally, the diffusive dissipation is stronger 
for a larger azimuthal gradient. The latter grows, on the average, linearly with N , that is 
why R min also increases with N , as seen in Fig. 2(a). 

If the initial radius of the necklace exceeds the above-mentioned maximum values, i.e., for 
the case of M = N, and ~ 1.8N for M ^ N, as per Eq. (3) (recall these two cases correspond 
to the fusion of the necklace into the fundamental and vortex soliton, respectively), Fig. 3 
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shows that the necklace arrays with M = do not undergo the fusion, but feature slow 
expansion, while their counterparts with M / make the number of the "beads" smaller 
via the fusion involving some of them, after passing several hundred diffraction lengths. The 
number of the lost "beads" is larger for smaller values of \M — N\, and smaller initial radius 
of the necklace, R . 

Finally, if the initial radius exceeds a still larger threshold value, i?^ n , the interaction 
between the "beads" in the necklace array becomes negligible, irrespective of its topological 
charge M (which implies R^ in is defined with some uncertainty, as it may slightly vary with 
the change of the total propagation distance). As a result, the pattern keeps its necklace-like 
structure and the initial radius, as shown in Fig. 4. For the same reasons as mentioned above 
[cf. Fig. 1(a)], the dependence of the respective minimum radius R^ in on modulation number 
iV of initial pattern (2) is approximately linear, see Fig. 4(a). Note that each individual 
element in the necklace patterns observed in Figs. 4 features an isotropic (circular) shape, 
unlike the "beads" in the initial pattern. This is explained by the fact that each element 
evolves into a fundamental soliton, whose deformation induced by the interaction with its 
neighbors is negligible, since the distance to them is too large. In fact, the "freezing" of 
soliton necklaces of a large radius due to the exponential decay of the interaction forces was 



observed in many cases before 
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IV. CONCLUSION 



We have demonstrated that, unlike conservative systems described by multidimensional 
nonlinear Schrodinger equations, models of dissipative optical media (in particular, laser 
cavities [3]), based on the 2D CGL equation with the cubic-quintic nonlinearity, admit 
fusion of necklace-ring patterns into stable fundamental and vortical solitons, provided that 
the initial radius of the ring, Rq, is not too large (if Ro is somewhat larger, the necklace 
may undergo a partial fusion). The outcome of the evolution is controlled by values of 
modulation number N and topological charge M of initial necklace (2). A simple analysis 
has predicted the vorticity of the fused soliton to be \N — M\, which is fully confirmed by 
direct simulations. Threshold characteristics of the fusion process, such as those displayed in 
Figs. 1(a) and 2 (a, b), were also explained in a qualitative form. The possibility to control 
the outcome of the fusion by means of the initial topological numbers, N and M, suggests a 
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principal possibility to use the process in all-optical switching schemes. On the other hand, 
if the initial radius of the necklace array is too large, the interaction between individual 
elements is negligible, allowing the pattern to keep its initial radius and necklace structure, 
while each element assumes the isotropic shape, corresponding to fundamental dissipative 
solitons in 2D. 

In addition to the straightforward realizations in terms of nonlinear optics, the model may 
also find application to Bose-Einstein condensates, where the dissipation and gain occur in 
models of matter-wave lasers, as demonstrated experimentally and theoretically in various 



settings [28, 



29 



30 



31 



32 



33 



34. 



35l . l36l . l37| . Thus far, such models were analyzed only 



in the effectively ID geometry, while the present results suggest to extend them into two 
dimensions. 
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FIG. 1: (Color online) The fusion into a fundamental soliton of the necklace array whose initial 
radius is not too large, Rq < R^ mx , and the topological charge is equal to half the number of 
"beads" in the array, M = N. (a): The largest radius, admitting the fusion, versus N. (b, c): 
Examples of the fusion for Rq = 4 and M = N = 5 (b) or M = N = 3 (c) [the examples are shown 
by means of contour plots of the local power, \u(x,y)\ 2 ]. 
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FIG. 2: (Color online) The fusion into a stable vortex soliton of the necklace whose initial radius 
and topological charge satisfy constraints R^m < Ro < 1.8iV, and either M^ in < M < N or 
N < M < M^ax — 2iV — M^ in . (a): The minimum initial radius admitting the fusion versus iV; 
(b): the minimum topological charge versus the initial radius for different fixed values of modulation 
number N. (c-e): Examples of the formation of the vortex soliton with N = 5 and Rq = 5 for 
M = 4 (c), M = 3 (d), and M = 2 (e). Additionally, panel (f) displays an example of the decay 
of the necklace cluster with N = 5 and M = 4, in the case of Rq < R^in (here, Rq = 1.9 and 
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FIG. 3: (Color online) (a): Slow expansion of the necklace array with N = 5, M = and i?o = H> 
in the case when its initial radius slightly exceeds the maximum value 1.8N, see Eq. (3). (b): 
In the same case, but with M = 3, ten initial "beads" fuse into eight and, eventually, into six 
individual elements. 
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FIG. 4: (Color online) Formation of "frozen" patterns which look like stable soliton necklaces, 
with the initial necklace radius exceeding R^ in and an arbitrary value of the topological charge 
(M). (a): iJ* n as a function of modulation number N (b, c): Examples of the formation of stable 
necklace rings with N = 5 and Rq = 12 for M = 2 (b) and M = (c). 
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